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vn : 

' We present the results of a search for high energy extraterrestrial neutrinos with the Baikal underwater Cherenkov 

(~| . detector NT200, based on data taken in 1998 - 2002 (1038 live days). Upper limits on the diffuse fluxes of 

Qh' i'e + t'/i + i^T, predicted by several models of AGN-like neutrino sources, are derived. For an E~'^ behavior of 

Q . the neutrino spectrum, our limit is E^<^^{E) < 8.1 x 10^^ cm^^ s^^ sr^^ GeV over an neutrino energy range 

$-H ' 2 X 10^ 5 X 10^ GeV covering 90% of expected events. The upper limit on the resonant v^. diffuse flux is 
to! <3.3 X 10-20 cm-2s-isr-iGeV~i. 



1. Introduction 

High energy neutrinos are likely produced in many violent processes in the Universe. Their detection would 
unambiguously reveal the hadronic nature of the underlying processes. Neutrinos would be generated by 
proton-proton or proton-photon interactions followed by production and decay of charged mesons. 

A description of the Baikal-detector as well as physics results from data collected in the years 1998 - 2000 have 
been presented elsewhere fT"2l|3]lll|5l. In this paper we present new results of a search for diffuse neutrinos 
with energies larger than 10 TeV. The analysis is based on data taken with the Baikal neutrino telescope NT200 
in the years 1998-2002. Instead of focusing to particles crossing the array, the analysis is tailored to signatures 
of isolated high-energy cascades in a large volume around the detector. This search strategy dramatically 
enhances the sensitivity of NT200 to diffuse high energy processes. 

The cascades can stem from leptons and hadrons produced in high energy charged current processes. Obvi- 
ously, the energy released by the hadronic cascade in NC-reactions is small compared to that of the leptonic 
cascades in CC-reactions. Since only electrons and taus develop cascades (the one by directly showering up, the 
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Figure 1. Left panel: the distribution of experimental events (dots) which survive condition Q as well as the expected 
distribution of simulated background events (histogram). Right panel: the N^^^ distribution of experimental events (dots) 
as well as background prediction (histogram). 

other via its decay to secondary particles which develop a cascade), the sensitivity of this search is dominated 
by Ve and ly-r detection. 

2. Data selection and analysis 

Within the 1038 days of the detector live time between April 1998 and February 2003, 3.45 x 10* events with 
^hit — 4 have been recorded. For this analysis we used 22597 events with hit channel multiplicity >15 
which obey the condition: 

*min = min(ti - tj) > -10 ns, i < j. (1) 

Here, t^, tj are the arrival times at channels i, j and the numbering of channels rises from top to bottom along 
the string. 

Figure^shows the t,T,in and N^^^ distributions for experiment (dots) and background simulation (histograms). 
The distribution of experimental events is consistent with the background simulation. No statistically signifi- 
cant excess above the background from atmospheric muons has been observed. 

With no experimental events outside the area populated by background events in the (t^i^, A'jjif)-parameter 
space, we derive upper limits on the fluxes of high energy neutrinos as predicted by different models of neutrino 
sources. 

The detection volume V^ff averaged over all neutrino arrival directions, rises from '^10^ m'^ for 10 TeV up to 
(4-6) X lO*' m^^ for 10'' TeV and significantly exceeds the geometrical volume Vg w 10^ m^^ of NT200. This is 
due to the low light scattering and the nearly not dispersed Ught fronts from Cherenkov waves originating far 
outside the geometrical volume. 

Since no event has been observed which fulfills the selection conditions, upper limits on the diffuse flux of 
extraterrestrial neutrinos are calculated. For a 90% confidence level an upper limit, ngo% =2.5, on the number 
of signal events is obtained according to Conrad et al. (6| with the unified Feldman-Cousins ordering 1 7 1 . 
We assume an uncertainty in signal detection of 24% and a background of zero events (which leads to a 
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conservative estimation of ngo% according to the Feldman-Cousins approach). If the expected numbers of 
signal events A'^^o^jgi is larger than ^90%, the model is ruled out at 90% CL. Table [^represents event rates and 
model rejection factors (MRF) n9o%/A^niodel for models of astrophysical neutrino sources obtained from our 
search. Recently, similar results have been presented by the AMANDA collaboration |8 9 1. Model rejection 
factors obtained by AMANDA are also shown in Table The models by Stecker and Salamon 1.1 OJ labeled 

Table 1. Event rates and model rejection factors for models of astrophysical neutrino sources. The assumed upper limit on 
the number of signal events with all uncertainties incorporated is ngo% = 2.5 



BAIKAL 


AMANDA 1 8 9| 


Model 










?l90%/^model 


?^90%/^model 


lO"*^ X E-'^ 


1.33 


0.63 


1.12 


3.08 


0.81 


0.86 


SS Quasar (To) 


4.16 


2.13 


3.71 


10.00 


0.25 


0.21 


SPu fTTl 


17.93 


7.82 


14.43 


40.18 


0.062 


0.054 


SPlfTTl 


3.14 


1.24 


2.37 


6.75 


0.37 


0.28 


P»^fT2l 


0.81 


0.53 


0.85 


2.19 


1.14 


1.99 




0.29 


0.22 


0.35 


0.86 


2.86 


1.19 


MPRitrti 


0.25 


0.14 


0.24 


0.63 


4.0 


4.41 


SeSi fT5l 


0.47 


0.26 


0.44 


1.18 


2.12 





"SS Q", as well as the models by Szabo and Protheroe 1 1 1 1 "SP u" and "SP 1" represent models for neutrino 
production in the central region of Active Galactic Nuclei and are ruled out with ?^9o% /^model ~ 0.06 - 0.4. 
Further shown are models for neutrino production in AGN jets: calculations by Protheroe 1121 . by Mannheim 
(HI, by Mannheim, Protheroe and Rachen 1 14 1 (model "MPR") and by Semikoz and Sigl iT^I "SeSi". The 
latter models for blazars are currently not excluded. For an E^"^ behaviour of the neutrino spectrum and a 
flavor ratio v f. : v : = \ : \ : \, the 90% C.L. upper limit on the neutrino flux of all flavors obtained with 
the Baikal neutrino telescope NT200 (1038 days) is: 

£;2$ < 8.1 X IQ-W^^s^isr^^GeV. (2) 

Assuming an upper limit on the number of signal events 71^0% =2.5, the model-independent limit on i7g at the 
W - resonance energy is: 

< 3.3 X 10"2"cm"2g-igj.-iQgY-i (3) 

Figure|2(left panel) shows our upper limit on the {v^ + v^ + Vr) E^^ diffuse flux as well as the model indepen- 
dent limit on the resonant flux (diamond). Also shown are the limits obtained by AMANDA and MACRO 
l|8| l9l 1161 . theoretical bounds obtained by Berezinsky (model independent (B) and for an E^"^ shape of the 
neutrino spectrum (B(i;"2)) fT7 1, by Waxman and Bahcafl (WB) fltl, by Mannheim et al.(MPR) 1 14 1, predic- 
tions for neutrino fluxes from topological defects (TD) 1 15 |, prediction on diffuse flux from AGNs according 
to Nellen et al. (NMB) 1191 . as well as the atmospheric conventional neutrino fluxes 1201 from horizontal and 
vertical directions (upper and lower curves, respectively) and atmospheric prompt neutrino fluxes obtained by 
Volkova et al. 1 2 1 j . Our upper limits (solid curves) on diffuse fluxes from AGNs shaped according to the model 
of Stecker and Salamon (SS) 1101 and of Semikoz and Sigl (SeSi) 1151 are shown in the right panel of fig. |2] 

In March/ April 2005 we fenced a large part of the search volume with three sparsely instrumented strings 
(see |4| for details). The three-year sensitivity of this enlarged detector A^r200+, with about 5 Mton enclosed 
volume, is approximately E^^^^ ^ lO^^cm^^s^^sr^^GeV for E > 10^ TeV, i.e. three-four times better than 
NT200. NT200+ will search for neutrinos from AGNs, GRBs and other extraterrestrial sources, neutrinos from 
cosmic ray interactions in the Galaxy as well as high energy atmospheric muons with i?^ > 10 TeV. 
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Figure 2. Left panel: neutrino flux predictions in different models of neutrino sources compared to experimental upper 
limits to E^^ fluxes obtained by various experiments, (see text). Right panel: experimental limits compared to two model 
predictions. Dotted curves: predictions from model SS 1101 and SeSi 1151 . Full curves: upper limits to spectra of the same 
shape. Model SS is excluded (MRF=0.21), model SeSi is not (MRF=2.12). 
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